Realistic many-body models for Manganese Monoxide under pressure 
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In materials like transition metals oxides where electronic Coulomb correlations impede a de- 
scription in terms of standard band-theories, the application of genuine many-body techniques is 
inevitable. Interfacing the realism of density-functional based methods with the virtues of Hubbard- 
like Hamiltonians, requires the joint ab mitio construction of transfer integrals and interaction ma- 
trix elements (like the Hubbard U) in a localized basis set. In this work, we employ the scheme 
of maximally localized Wannier functions and the constrained random phase approximation to cre- 
ate effective low-energy models for Manganese monoxide, and track their evolution under external 
pressure. We find that in the low pressure antiferromagnetic phase, the compression results in an 
increase of the bare Coulomb interaction for specific orbitals. As we rationalized in recent model 
considerations [P/ii/s.i?e?;._B79, 235133(2009)], this seemingly counter-intuitive behavior is a conse- 
quence of the delocalization of the respective Wannier functions. The change of screening processes 
does not alter this tendency, and thus, the screened on-site component of the interaction - the Hub- 
bard U of the effective low-energy system - increases with pressure as well. The orbital anisotropy 
of the effects originates from the orientation of the orbitals vis-a-vis the deformation of the unit-cell. 
Within the high pressure paramagnetic phase, on the other hand, we find the significant increase of 
the Hubbard U is insensitive to the orbital orientation and almost exclusively owing to a substantial 
weakening of screening channels upon compression. 



PACS numbers: 

I. INTRODUCTION 
A. Realistic models for correlated materials 

While standard band structure methods, like density 
functional theory (DFT) in the local spin density ap- 
proximation (LSDA),^'^ in principle allow for a treatment 
of the realistic complexity of materials arising from or- 
bital, spin and structural degrees of freedom, the failure 
of these techniques in the presence of strong electronic 
correlations is well documented.'^ Many-body techniques, 
on the other hand, while potentially treating correlation 
effects with accuracy, often lack the realism to account 
for the electronic structure of a material from a first- 
principles basis. 

Hence, when it comes to the ab initio description of 
correlated matter, the joining of these two fields is of 
paramount importance. Often, this combining amounts 
to constructing effective, yet realistic, low energy models, 
in the spirit of the Hubbard or Anderson model. These 
consist of a one-particle part, extracted from the Kohn- 
Sham band structure, and two-particle interaction terms. 
There are different paradigms for choosing the matrix el- 
ements of the latter. Sometimes, these values are treated 
as mere empirical parameters, adjusted to correctly re- 
produce some experimental finding. However, there exist 
first-principle methods, among which the most popular 
are the constrained LDA technique,^ and the constrained 
random phase approximation (cRPA).^ 

In this article we are chiefly concerned with the evolu- 
tion of interaction matrix elements under pressure within 
a realistic setup. While the influence of pressure onto 



structures and band-structures has been widely studied 
on the DFT level, changes in the interaction have re- 
ceived much less attention. Yet, the pronounced sensi- 
tivity of correlated matter on external parameters^ her- 
alds the need for an accurate treatment of all ingredients 
to the electronic structure. For the LDA-|-U^ technique, 
this issue has been investigated within a linear response 
based approach*"^", and the reproduction of crystal and 
electronic structures, as well as phase stabilities were in- 
deed found to require the description of pressure induced 
changes in the Hubbard U . 

Here, we are interested in setting the stage for many- 
body methods that go beyond density functional meth- 
ods. In this vein, a particularly versatile technique is the 
use of maximally localized Wannier functions^^'^^ within 
the cRPA setup, ^"^'^^ which is the method of choice of 
this article. In this formalism it is possible to turn off 
precisely those screening channels that are to be left to 
the solution of the effective model. ^'^^ Moreover, the use 
of many-body techniques often requires the use of local- 
ized basis sets, and the Wannier construction is geared 
at exactly this. While the Wannier-cRPA technique has 
been already applied to several systems, such as 3d tran- 
sition metals, the perovskite SrVOa,^^ and oxypnic- 
tides"'^^^^^, the evolution of the respective Coulomb in- 
teraction matrix elements under external stimuli such as 
pressure have not been addressed in the realistic case. 
In recent model considerations,^^ we established generic 
behaviors of maximally localized Wannier functions and 
interaction matrix elements under pressure for simple se- 
tups. We rationalized that as a consequence of the delo- 
calization of Wannier functions under compression, not 
only do transfer integrals increase, but also - counter- 
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intuitively - the matrix elements of the bare Coulomb 
interaction augment. In this paper, we extent this work 
to a more complex scenario by applying the fully ab ini- 
tio Wannier-cRPA formalism to manganese monoxide, 
MnO, under pressure. In particular, we monitor the 
evolution of transfer integrals, on-site matrix elements 
of the bare Coulomb interaction, and of the partially 
screened Coulomb interaction (the Hubbard U, and the 
Hund J), for different realistic low-energy systems. We 
will distinguish between pressure induced changes related 
to the construction of the sub-system, and changes in 
the screening strength as provided by the band-structure 
and transition matrix elements in the polarization. We 
find that while the former mainly determine the behav- 
ior of the interactions in the low temperature, low pres- 
sure anti- ferromagnetic phase of MnO, changes in the 
Hubbard U of the high pressure paramagnetic phase are 
dominated by screening related effects. 



B. Manganese monoxide 

1. Approaches to tfie crystal and electronic structure 

Concerning the failures of the LDA in the realm of cor- 
related materials, the focus material of this work, man- 
ganese monoxide, MnO, is rather well behaved at first 
glance : Below 118 K the system is antiferromagnetic 
and due to the fact that the Mn3d orbitals are half- 
filled, spin-dependent band-structure calculations ^^"^"^ do 
indeed yield an insulator. However, above the Ncel tem- 
perature, which is much smaller than the experimental 
charge gap of 3.9 eV,^^ MnO is paramagnetic, and be- 
lieved to be a Mott insulator. A straight forward band- 
structure calculation in the paramagnetic phase, how- 
ever, yields a metal. 

A genuine feature of correlated materials is their sensi- 
tivity to external stimuli, leading to rich phase-diagrams 
with regimes of different spin or orbital orders, and metal- 
insulator transitions.^ MnO is no exception and both ex- 
periments and theory helped to elucidate the phases of 
manganese monoxide as a function of temperature versus 
pressure. See Fig. 1 in Rcf.^^ for a graphical depiction 
of compiled results. 

The paramagnetic insulator MnO, which is of rock- 
salt structure (Bl), contracts below the Neel tempera- 
ture along the (111) direction, resulting in a rhombohe- 
dral distorted phase (dBl). MnO has been the prototype 
material for the discovery and description of the superex- 
change mechanism. ^''"^^ It was rationalized by Terakura 
et al}'' that this strong manganese superexchange via the 
oxygen orbitals^'"'^® is responsible for a spin oder that is 
of AF II type,^""'^^ with alternating spins along the dis- 
tortion axis. With increasing pressure, the Neel temper- 
ature rises^'^''^^ and reaches room temperature (300 K) at 
30GPa.^^ It has been shown that beyond the volume com- 
pression, it is the rhombohedral distortion that stabilizes 
the AF II order under pressure with respect to other spin 



orientations. -'^^ At constant temperature (300 K), a fur- 
ther transition occurs at 90 GPa,^^'^^''^^ when the system 
becomes, again, a paramagnetic insulator, yet of NiAs 
(B8) structure. Notably, the high pressure B8 struc- 
ture was predicted by DFT calculations in the generalized 
gradient approximation (GGA).^^ At still higher pres- 
sure, the system undergoes an iso-structural Mott insu- 
lator to metal transition at 105 GPa,^^ signaled by a vol- 
ume collapse of 6.6%,^^ metallic lustre, "^^ an increase of 
reflectance, ^^'^^ a resistance decrease of several orders of 
magnitude^^ and a quenching of the magnetic moment. 
The transition was anticipated, from theoretical consid- 
erations'"^ and ah initio calculations.^^ Yet, both works 
forward different mechanism : In Ref.^° a pressure in- 
duced increase in the crystal-field was invoked to explain 
the mctalization, while in Ref.^^, in which the rhombo- 
hedral distortion is neglected, a growing band-width was 
seen at the origin of the transition. We will come back 
to this in the discussion of our results, which favor the 
first point of view. 

On the theory side, numerous works that aim at 
reproducing the described phase-diagram from first 
principles resort to improvements within the den- 
sity functional formalism. Among these are the 
generalized gradient approximation (GGA),^*"^"' self- 
consistent Hartree,*! LDA/GGA-|-U,i8'20-23,42,42 ggjf, 

interaction corrected (SIC) LSDA.'*'^''*^ Also applied were 
various hybrid functional"*^ based approaches^'^'^*'^'^'^''^. 

Although relying on the coincidence of a half-filled d- 
shell and the spin order for insulating results, it is due 
to the fact that (unlike for the other transition metal 
monoxides) the 3d electrons are not falsely itinerant, that 
the above methods succeed in yielding at times reason- 
able structures, lattice constants, bulk moduli, magnetic 
moments, and sometimes even values for the charge gap 
(see especially the SIC and hybrid functional references 
from above). Results have been obtained also within the 
GW approximation.^^*® While pioneering work*^ still 
had to resort to numerical simplifications, the increase in 
computer performance nowadays even allows the applica- 
tion of self-consistency schemes to MnO,^° which consid- 
erably improve on the gap value of the antiferromagnetic 
phase with respect to LDA calculations. 

Yet, despite these successes, it is evidently desirable 
to have a many-body technique that works in a wider 
context, capturing correlated metals, ordered and Mott 
insulators. In this vein, the probably most generally ap- 
plicable approach nowadays is LDA-|-DMFT,^* where a 
Hubbard Hamiltonian consisting of a one-particle part 
stemming from the LDA and an interaction part of 
Hubbard-Hund type is solved within dynamical mean 
field theory (DMFT).52,68 ^j^^g Qprp j^ethods and 

their generalizations, including LDA-|-U,^ cannot ad- 
dress finite-temperature insulators without broken sym- 
metries, LDA-I-DMFT is in principle capable to deal with 
all the different phases. Concerning MnO, the tech- 
nique has first been used for the calculation of phonon 
dispersions.^'^ Recently LDA+DMFT was applied by 
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Kunes et al. to study the isostructural Mott insula- 
tor to metal transition witnessed at lOSGPa (300 K). 
The results reproduced quite accurately all experimental 
findings, in particular the critical pressure and the joint 
collapse of volume and magnetic moment. 

Often it is criticized that LDA+t/ and LDA+DMFT 
calculations are not fully ah initio since the interactions 
are treated as "parameters" that arc in a sense tuneable. 
Indeed in the above LDA+DMFT work,^'' the Hubbard 
U and the Hund J were chosen to be constant as a func- 
tion of pressure. Although adjustable parameters are to 
be preferred over artificially broken phases, it is true that 
especially in transition metal monoxides there are several 
competing energy scales^^ so that the values of the inter- 
actions have to be known with greater precision than in 
other compounds. Among those relevant energies are : 
the crystal field splitting, that measures the degree of 
non-degeneracy of the centers of gravity of the d orbitals, 
e.g. Ac/ = Egg — et2g for 8-11 octahedral coordination, the 
charge transfer energy Act = ^t2g ~ ^2p here defined be- 
tween the 02p and the Mn t2g, the transfer integrals 
(hoppings) tf^g, tcg , the hybridizations tpd, and the on- 
site Coulomb interaction U of the d orbitals. 



2. Earlier works on the Hubbard U 

In the pioneering work of Anisimov et alJ that intro- 
duced the LDA-I-C7 formalism, the constrained LDA tech- 
nique was applied to MnO at zero pressure and a value 
U = 6.9 cV was found. The Hubbard U is thus of the 
same order as the charge transfer energy Act (see below) . 
Therewith MnO is in the intermediate regime between 
a charge transfer and a Mott-Hubbard insulator^^ and 
changes in U might considerably alter the many-body 
spectra and their interpretation. 

An empirical technique that was applied to MnO is the 
fitting of configuration interaction (CI) cluster calcula- 
tions to photoemission experiments. Although of course 
depending on the details of the chosen cluster, we cite the 
values U = 7.5 eV, Act = 7.0 eV of Fujimori et al.^*^ and 
U = 8.5 eV, Act = 8.8 eV of van Elp et al..'^^ Still, CI 
cluster calculations cannot account for delocalized bands. 
As a consequence, the struggle of d electrons between lo- 
calized and itinerant behavior and the mixing with the 
more extended 02p and Mn4s orbtials is biased. 

Using a semi-empirical Anderson (d-) impurity Hamil- 
tonian, Zaanen and Sawatzky^^ estimated A and U to 
be both around 9 eV, and found accurate values also for 
the Ncel temperature. An Anderson Hamiltonian was 
recently also applied to results from resonant inelastic 
x-ray scattering (RIXS) spectra^^ : the parameters of a 
single manganese d impurity model were adjusted to yield 
U ~ 7.2 eV and Act = 6.5 eV. In this approach, the de- 
localized character of the oxygen and Mn4s orbitals is 
better captured than within CI based methods. Indeed, 
this is a step toward LDA+DMFT calculations, albeit 
without obeying the self-consistency condition that re- 



lates the manganese impurity with the periodic solid. 

The Hubbard U can also be extracted from GW calcu- 
lations. Indeed the RPA technique is a crucial ingredient 
to the GW approximation. Pioneering model GW calcu- 
lations suggested a value oiU = 8.0 eV*^ for the antifcr- 
romagnetic phase, albeit only empirically extracted from 
the energy difference between occupied and unoccupied 
Cg and t2g states. 



II. METHOD 

The method of using the Wannier orbital construction 
in conjunction with the constrained RPA has been pre- 
sented in detail in Ref.^^. Here, we shall just give a brief 
summary as is necessary for understanding our results. 



A. Band-structure 

In this work, we employ the LSDA in the full-potential 
(FP) LMTO^^ realization and use experimental struc- 
tures and lattice parameters as provided by Yoo et al?^ 
at 300 K®^ and variable pressure. We use up to 4096 
reducible Brillouin zone points, the LMTO basis set in- 
cludes orbitals up to ^ = 4(3) for Mn(0), and we use local 
orbitals^* to incorporate semi-core and high-lying states. 



B. Maximally localized Wannier functions 

The one-particle part of the effective low energy sys- 
tems is extracted from the (LMTO) band-structure cal- 
culation by the construction of a Wannier basis^^ of 
the desired subset of orbitals. This procedure not be- 
ing unique, we choose to use Wannier functions that 
are maximally localized. ^^'^^ While the initial Hilbert 
space is spanned by all Kohn-Sham wave functions, 
% — span{?/'^^}, for the aim of constructing effective 
low energy models necessitates the choosing of a sub- 
space, 'U'-°^ c %, made up by only selective Kohn-Sham 
orbitals. The Wannier functions are then deduced from 
the latter. Complications occur if the bands of the or- 
bitals that constitute the desired sub Hilbert space are 
entangled with high energy bands. For details on how 
then to construct maximally localized Wannier functions 
see^^, and for a comparison with other Wannier function 
schemes, for instants the approach in which the on-site 
screened Coulomb interaction is maximized, see e.g. the 

Rcfs.13.14.60, 



C. constrained RPA 

The matrix elements of the bare Coulomb interaction 
f (r, r') = 4fj^l/ |r — r'| in the Wannier basis XRa('") ^'"^ 
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setup 


low energy subspace rt 


polarization Piow 


d 

pd 
pd+p 


span {i)ka 01 = Mn3d} 
span {^^^ 1 a = Mn3d, 02p} 
spanlV;^/ 1 a = Mn3d, 02p} 


Pd 

Pd+Pp+Ppd+Pdp 
Pd 



TABLE I: Summary of the low-energy setups that are used : 
Shown are the choices for the construction of the low energy 
subspace, and the polarization. See text for details. 
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(1) 



From the latter, partially screened interaction matrix el- 
ements are obtained using the cRPA technique. ^■^'^ Using 
the Kohn-Sham orbitals, the total RPA polarization can 
be expressed by 



P(r,r',L.) = 



(2) 



occ unocc 



EE E V'J^r(r)^K.'(r)«(r')C;; (r') 

spin kn k'ri' 

1 1 



W - Ck'ri' + fkn + lO+ UJ + ek'n' - Ckn - «0+ 

Its matrix elements within the Wannier basis set are 
given analogous to the above bare interaction. The 
strength of screening channels can thus be influenced by 
the overlap integrals (matrix elements) of the wave func- 
tions, and by the energy difference of the Kohn-Sham 
excitations that appear in the denominators. 

From the form of the fully screened interaction, 
W{u!) = [1 — P{uj)v]^^ V (that also appears in the GW 
formalism"*^^^* ) it can be shown, ^ that the screening con- 
tributions are actually additive. Indeed one can split the 
polarization, P = Plow + Pr, into the transitions within 
a low-energy orbital subspace. Plow of and the 

rest Pr = Phigh-high + Phigh-low + Plow-high, wherc the 
latter includes, both, transitions within the high-energy 
subspace, and between the two subspaces. Then the 
fully screened interaction can be written recursively, i.e. 
W ^[l- PiowWl""^] ~ V/°"' where 



1 - PrV 



(3) 



is the Coulomb interaction when all screening processes 
that do not involve transitions within the low-energy or- 
bital subspace have been accounted for^" . 

For our compound, MnO, we will consider three differ- 
ent setups, that are also summarized in Tab. I. 



• In the "d-setup" , a many body model is obtained 
from maximally localized Wannier functions that 
are constructed for the space of the MnSd Kohn- 

and from a partially 



Sham (KS) orbitals, V^f^, 



screened interaction Wf that is given by constrain- 
ing the polarization to transitions that are not fully 
within the d-subspace^^'^^. 

• MnO being close to the charge transfer regime, it 
is sensible to include the oxygen 2p orbitals for 
a realistic many-body treatment. In such a "pd- 
setup" one constructs a Wannier basis for the space 
of MnSd and 02p Kohn-Sham wavefunctions, and 
chooses Pr ^ P - Pd~ Pp- Ppd- Pdp- 

• Yet, often, a Hubbard U is only put on the most 
correlated orbitals, here the Mn3d ones. In that 
case, while still constructing a full pd Wannier ba- 
sis, one should allow for oxygen screening in the 
RPA (i.e. Pr = P — Pd), since those channels are 
blocked when the interactions between the p and 
d orbitals are omitted. Therewith, Wr(r, r') of the 
current setup and the pd one are the same, yet 
their matrix-clement will differ due to the Wannier 
different basis. 

As an outlook, we like to mention that in the described 
procedure to generate the effective many-body problem, 
couplings in the interaction between the low energy sub- 
space and the other orbital degrees of freedom arc dis- 
carded - an approximation that might prove insufficient 
if the subspace is chosen too restrictively. A framework 
that separates the Hilbert space, while retaining the in- 
fluence of self-energy effects between subspaces has been 
recently proposed in Ref.^^. 



D. Wannier functions and Coulomb matrix 
elements under pressure 

Before applying the above techniques to our compound 
of interest, we shall first discuss what circumstances will 
influence the Wannier functions and the interaction val- 
ues. We can distinguish three conceptually different 
mechanisms : The choice of the orbital sub-space, the 
consequences of structural changes onto the localization 
of Wannier functions, and the changes related to the 
strength of screening processes. 



1. The choice of the orbital sub-space 

The construction of an effective low energy system de- 
pends very much upon the choice of the orbitals that it 
consists of. In particular, the dimension of the orbital 
sub-space 'H'""' has a strong impact on the localization 
of the Wannier functions. Indeed, from the maximally lo- 
calized Wannier function point of view, a larger subspace 
will allow, for a given orbital, for the mixing in of more 
Kohn-Sham orbitals. Therewith, the variational freedom 
to localize the basis functions becomes larger, and the 
interaction elements bigger. See, for an example, Ref.^^. 
We thus expect the interaction matrix elements of the 
pd-model to be larger than in the d-only model. 
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2,. Wannier functions under pressure 



External pressure, or structural changes in general, 
provide another impetus to alter not only the one-particle 
band-structure, but also the Wannier functions and thus 
the bare Coulomb interaction matrix elements. 

The fact that Wannier functions are not eigenfunc- 
tions even of the non-interacting problem, may result in 
counter- intuitive tendencies. In Ref.^^, we established 
generic behaviors of maximally localized Wannier func- 
tions and interaction matrix elements under pressure. In 
order to put things into perspective, we shall here briefly 
summarize some findings of this work : 

As a model system, we consider a tight-binding 
parametrization of a one dimensional solid, in which 
case Wannier functions are already maximally localized 
if inversion symmetry is verified, ^^'^^ and limit the dis- 
cussion to the case of a single band. Details can be 
found in Ref.^**. As building blocks of the basis func- 
tions we take hydrogcn-like Is orbitals in one dimen- 
sion, x(x) = l/y^oo exp(— |x| /oo); with the Bohr ra- 
dius ao. The Bloch eigenfunction is given by tpk{x) = 
AkY.R(^'''^x{x - R), where Ak assures the normaliza- 
tion. The maximally localized Wannier functions are 
then iPb{x) = "^n-^nXi.^ — na — R), with the lattice 
constant a, while An is a Fourier transform of Ak and 
also incorporates the normalization of the Wannier func- 
tion. In the limit of large atomic separation (a 3> ao), the 
Wannier functions will equal the atomic orbitals. Thus 
An = Sn,o for the representative site "0". When pres- 
sure is applied, and the lattice constant shrinks, overlaps 
of the (non-orthogonal) atomic orbitals causes contribu- 
tions from neighboring sites to mix into the basis func- 
tions, and the distribution An broadens. This leads to an 
increase in the (nearest neighbor) transfer integral t and 
a concomitant increase of the spread (as defined in^^). 
In this sense, charge carriers become more dclocalized, 
as expected. The tails of the Wannier functions extent 
over several lattice constants, before the originally expo- 
nential decay sets in. However, as can be shown, the 
coefficient Aq in the above decomposition on the array 
of atomic orbitals, i.e. the strength of the on-site orbital, 
grows, and thus becomes larger than one. This results 
in a greater probability density |'0_r(x)| around the site 
origins, x — R — 0, when pressure is applied. This in turn 
causes a larger local Coulomb interaction matrix element 
V, Eq. (2) . Hence, the pressure induced delocalization is 
accompanied by a larger local interaction, whereas, intu- 
itively, one might have expected the opposite. In our one 
band toy-model, since Pr = (see Eq. (3)) the screened 
Coulomb interaction U equals the bare Coulomb interac- 
tion V. Moreover it is maximal in the maximally local- 
ized Wannier basis. 



P [GPa] II tp-t,^ I Act I / Aef I tt,, ^t,, 

(dBl)|| 0.60 I 3.5 I OlO | 0.21 
100 (dBl) 0.83 5.5 0.13 0.20 
100 (Bl) 0.91 6.5 0.13 0.41 

TABLE II: Comparison between the direct and oxygen me- 
diated contributions to the effective transfer element of the 
occupied t2g orbitals for the Bl/dBl phase of MnO. Act = 
e?2g ~ ^2p is the charge transfer energy. Shown are the largest 
transfers for nearest neighbors of same spin. All energies in 
eV. 

3. Screening strength and the interaction. 

While the above changes in overlaps and hybridizations 
directly enter the calculation by modifying the Wan- 
nier functions, that influence the matrix elements of any 
quantity, pressure moreover enters the construction of 
the effective low energy system by means of the band- 
structure. Indeed, according to Eq. (3), pressure induced 
modifications in the one-particle excitations will alter the 
screening strengths of possible transitions. When band- 
widths, crystal fields or bonding/antibonding splittings 
get enhanced upon the compression of the system, the 
screening will become less effective and, therewith, the 
Hubbard U larger. Besides this energetic efi'ect, evident 
from the denominators in Eq. (3), there is also an effect of 
matrix elements when expressing the polarization in the 
Wannier basis. For an example of the interplay of both 
effects see Ref.^^. Contrary to the preceding section, the 
transition energies and the matrix elements affect not the 
bare Coulomb interaction, but only screened quantities. 



III. RESULTS AND DISCUSSION 

A. MnO in the low pressure regime : The 
antiferromagnetic insulating dBl phase 

1 . bandstructure 

Fig. 1 displays the LSDA band-structure of antifer- 
romagnetic MnO for different pressures and structures. 
Explicitly marked are also the bands constructed from 
the Wannier orbitals that are used for the construction 
of the effective "pd" low energy system (the bands of the 
effective "d-only" model achieve the same precision) . Us- 
ing the maximally localized Wannier function formalism 
for entangled bands, the choice of Hilbert space for the 
Mn3d and 02p orbitals results in a clean distinction from 
the Mn 4s band that crosses the d-bands. 

The corresponding LSDA densities of states (DOS) of 
the pd-model are shown in Fig. 2. As we can see, pres- 
sure induces an increase in the crystal-field splitting, Ac/, 
i.e. a larger separation between the Cg and Kohn- 
Sham excitations (both the occupied and unoccupied), 
to an extent that the moving of the Cg weight towards 



(a) OGPa - Bl (NaCl) (b) OGPa - dBl (c) 52.5GPa - dBl (d) lOOGPa - dBl (c) lOOGPa - Bl (NaCl) 



FIG. 1: (color online) LDA band-structure of antiferromagnetic Bl/dBl MnO. Shown are the Kohn-Sham bands (solid lines) 
and the eigenvalues of the Wannier Hamiltonian (dashed). The vertical axis is energy (in eV) with the Fermi level at the origin. 
While (b) - (d) correspond to the true distorted dBl structure at different pressures, (a) & (d) show the band-structure for the 
artificial Bl (NaCl) structure for the indicated pressures. 
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FIG. 2: (Color online) Density of states (DOS) (in 1/eV) of 
antiferromagnetic Bl/dBl MnO for increasing pressure (top 
to bottom). Shown are the site-summed contributions of one 
spin component of the pd-model. The lowest panel displays 
fictitious Bl MnO at lOOGPa at the experimental volume. 
Different partial DOS are depicted as indicated in the graph 
at the bottom. For P = OGPa, the DOS of the dBl and the 
Bl phase do not differ significantly on the shown scale. 



the Fermi level eventually causes the system to become 
metallic.*^ Concomitant with the increase in crystal- 
field, pressure results in a growing Mn-0 hybridization 
tpa = (pO 1 77"^^ I do), where 77"^^ is the Kohn-Sham Hamil- 
tonian. This effect is especially large for the Cg orbitals, 
as is clear from the octahedral coordination. For values. 



see the later discussed Tab. II (occupied t2g orbitals), and 
III {Cg orbitals). Most notably is the fact, that the partial 
d-bandwidths do not increase with pressure. The origin 
of this is identified when comparing results for the exper- 
imental dBl phase with the undistorted Bl (NaCl) struc- 
ture of the same volume. As evident from Fig. 1 (d,e), 
and Fig. 2, the omission of the rhombohedral distortion 
causes the bandwidths to substantially increase. Indeed, 
in fictitious Bl MnO, the closing of the gap can be mainly 
attributed to the change in bandwidth, as was stated in 
j^gf 19,63 -^yj^gj-e distortion was neglected. The evolu- 
tion of the bandwidths can be quantified by looking at 
the effective MnSd nearest neighbor transfer elements for 
a given spin. The latter consists of two major contribu- 
tions : the direct hopping tdd = {dO\H'^^\dIl) (R nearest 
neighbor unit-cell of 0), and the oxygen mediated trans- 
fer tp^/Act, with the hybridization tpd, and the charge 
transfer energy Ac*. The (largest) nearest (spin) neigh- 
bor transfers of the occupied t2g orbitals, tt.^ _t.-^ , the 
corresponding hybridization with the oxygen 2p orbitals, 
tp-t2g , and the charge transfer energy Act = ^t2g ~ ^2p a-s 
measured by the band centers^^, are collected in Tab. II 
for ambient pressure and lOOGPa, for both, the experi- 
mental dBl and the undistorted Bl structure. 

As a result of the joint increase oitpd and Act, the oxy- 
gen mediated contribution to the t2g bandwidth remains 
roughly of the same magnitude, irrespective of pressure 
and distortion. The direct t2g-transfer for the distorted 
dBl structure is also constant with pressure. However, 
it is substantially structure dependent. Indeed the near- 
est neighbor t2g-t2g transfer doubles when the distortion 
is neglected, causing the effective hopping to augment 
by 64%, demonstrating the critical need to include the 
distortion for a proper description of dBl MnO under 
pressure. As to the unoccupied t2g orbitals, the transfers 
and energy differences are larger, but the same argument 
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FIG. 3: (Color online) local Coulomb interaction of the antiferromagnetic dBI phase of MnO for different pressures, (a) 
diagonal elements of the local, bare interaction V for the pd and d-only model, and resolved for the different orbitals a. (b) 
zero frequency limit of the RPA partially screened local interaction Uaa = M^ro o'""(i^ = 0) for the pd model (with and without 
p screening), and the d-only model. For comparison are shown also the results of undistorted, Bl (NaCl) structured MnO. 



as above holds, for the trends are the same. Wc conclude 
that within LDA (in contrast to previous works^^ (see 
also the recent^^) relying on the undistorted structure), 
the gap closure under pressure is propelled by the change 
in crystal field and not by a broadening of band- widths. 



P [GPa] 


,unocc{occ) 


A unocc{occ) 




,unocc{occ) 


(dBI) 


1.24 (1.08) 


8.0 (4.75) 


0.19 (0.25) 


0.05 (0.06) 


100 (dBI) 


1.72 (1.50) 


10.5 (7.50) 


0.28 (0.30) 


0.13 (0.11) 


100 (Bl) 


1.82 (1.63) 


10.5 (7.50) 


0.32 (0.35) 


0.14 (0.13) 



TABLE III: Comparison between the direct and oxygen me- 
diated contributions to the effective transfer elements of the 
unoccupied (occupied) eg orbitals for the Bl/dBl phase of 
Mn. Act = Eeg — e2p is the charge transfer energy. Shown are 
the largest transfers for nearest neighbors with the same spin. 
All energies in eV. 



The relevant quantities for the unoccupied (occupied) 
Cg orbitals are compiled in Table III. For the occupied 
eg orbitals the situation is opposite to the t2g ■ While 
the oxygen mediated contribution to the band-width, 
tp_g_^/A°j^, at lOOGPa depends sizeably on the struc- 
ture (Bl/dBl), the direct hopping is rather insensitive 
in this respect. For the unoccupied Cg none of the val- 
ues at lOOGPa depends much on the structure, as can 
be inferred also from the similar shapes of the density of 
states. Yet, overall the effective bandwidths of the dBI 
phase are always a little smaller than for the undistorted 
phase. Thus, the rhombohcdral distortion is effectively 
reducing all partial d-bandwidths. 



2. The Coulomb interaction 

a. The bare interaction. As explained in a model 
context in Ref.^^, changes in the orbital overlaps and hy- 
bridizations are expected to have an impact on the bare 
Coulomb interaction matrix V of Eq. (2). The on-site el- 
ements of the latter for antiferromagnetic dBI MnO are 
shown in Fig. 3(a), both for a pd and a d-only Wannier 
basis set, and for different pressures. 

Owing to the distortion axis, the degeneracy proper- 
ties of the Mn3d orbitals do not change with pressure. 
Moreover, in the current case, the maximally localized 
Wannier functions retain the degeneracy of the LMTOs, 
and while the LDA wavefunctions were merely the initial 
guesses for the maximally localized Wannier functions, 
we keep the labeling of "t2g", "eg"j "rf^y" and alike. In 
some cases, especially for entangled bands, the localiza- 
tion procedure might cause the resulting Wannier func- 
tions to acquire a different symmetry. 

In the case of the pd-model, the matrix elements are 
sensitive to pressure only for the 02p and the unoccu- 
pied Cg orbitals. There, the changes in the interaction are 
quite significant and reach up to 10%. This is a natural 
consequence of our above remarks : It is precisely those 
orbitals whose bondings changes the most under pres- 
sure, namely the hybridizations between them increase 
considerably. As summarized in the Method section, and 
discussed in detail in ReL^'*, this allows the respective 
Wannier functions to accumulate more weight at small 
distance, while an overall increase in spread accounts for 
the gain in hybridization. This is also the reason why, 
for the pd-model, the eg elements are larger than the t2g 
ones : the larger hybridizations of the eg orbitals with 
the oxygen 2p ones, allow the corresponding maximally 
localized Wannier functions to reduce their extension. 

As anticipated, the elements of the d-model are lower 
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P [GPa] 


C/pd [eV] 


UZ [eV] 


0.0 (dBl) 


2.25 


1.51 


52.5 (dBl) 


2.44 


1.69 


100.0 (dBl) 


2.65 


1.96 


100.0 (Bl) 


2.61 


1.74 


96.5 (B8) 


1.66 


1.11 


120.0 (B8) 


1.86 


1.26 


160.0 (B8) 


1.98 


1.42 



TABLE IV: Averaged nearest neighbor inter-atomic interac- 
tions f/pd = Wrl^/^u = 0) and = Wri^^^'^^^^uJ = 0) 
(where di/2 is a short hand for the d-orbitals of the 1./2. Mn 
atom in the cell) for the pd-model of the Bl/dBl phase (top) 
and the B8 phase (bottom) for different pressures. 



than for the pd case, since the Wannier functions are 
more extended. Having restrained the variational free- 
dom to the d orbitals, their matrix elements acquire a 
larger dependence on pressure than in the pd case. This 
is a result of the fact that the intra-d hoppings, especially 
the Cg ones, must make up for the dispersion caused by 
oxygen mediated transfers that are not present in the 
smaller basis. 

h. The partially screened interaction. While the 
bare Coulomb interaction was mainly a mean to analyze 
the pressure dependence of the Wannier functions, the 
quantity that is needed for the formulation of the effec- 
tive models is the partially screened Coulomb interaction. 
In Fig. 3 (b) we show, for the same pressures as before, 
the Hubbard U for the three different types of models, as 



introduced in the Method's section : the full pd-model, 
the pd-model with pd-screening (for a treatment in which 
only the d-orbitals will be supplemented by a local inter- 
action), and the d-only model. The magnitudes of the 
corresponding interactions decrease in the given order : 
The interactions in the pd-model are smaller when the 
p-screening is present, and the interactions in the d-only 
model are yet smaller because of the larger extension of 
the Wannier functions. The pressure dependence of the 
elements has the same tendencies as the bare Coulomb 
interactions. Yet, we note that the relative difference, 
^jjWOGPa _ jjOGPayjjOGPa^ jg ^^^^^^^ fo^, ^j^g Hubbard U 

than for V, the pressure influence thus stronger. This is 
understood from a decrease in screening strength owing 
to the larger separation in energy between occupied and 
unoccupied states as pressure increases (see e.g. Ref.^^ 
for a discussion), as apparent from the band-structure, 
Fig. 2. On an absolute scale, and for the given pressure 
range, the Hubbard U changes by up to 1 eV for both 
the Mn3c? and the 02p orbitals ~ a non-negligible effect 
given the aforementioned modifications in the magnitude 
of the transfer integrals. 

For the sake of completeness and to put the values 
further into perspective. Table IV (upper part) con- 
tains the nearest neighbor partially screened interactions 



0) and 



nn 
dd 



= 0) 



(where di/2 is a short hand for the d-orbitals of the 1./2. 
Mn atom in the cell) within the pd-model and averaged 
over the respective orbital subset. Both elements are 
roughly a factor of four smaller than the intra-atomic 
interactions. 



(100) GPa 


dxy 


dxz 


d^2 


dyz 


dx"^ — y'^ 


dxy 




0.61 (0.62) 


0.66 (0.68) 


0.61 (0.62) 


0.38 (0.37) 


dxz 


0.61 (0.62) 




0.46 (0.45) 


0.61 (0.62) 


0.59 (0.60) 


d^2 


0.66 (0.68) 


0.46 (0.45) 




0.45 (0.45) 


0.58 (0.65) 


dy^ 


0.61 (0.62) 


0.61 (0.62) 


0.45 (0.45) 




0.59 (0.60) 




0.38 (0.37) 


0.59 (0.60) 


0.58 (0.65) 


0.59 (0.60) 





TABLE V: Hund's coupling J = Wro o (^ = 0,a 7^ /3) of antiferromagnetic dBl MnO and the d-only model for the extreme 
pressures and lOOGPa. 



As a function of pressure, both Upd and V^^ increase. As pressure), it does not depend significantly on pressure, 
was true for the diagonal elements of the Hubbard U, the 
pressure dependence of these off-diagonal interactions is 
a feature mostly intrinsic to the Wannier functions and 
not the screening strength. 

The Hund's coupling, J = W^,.o,o'^"(w = 0,a 7^ /3), 
is compiled below in Table V. While being largely 
anisotropic with respect to the orbital (for the Mn d or- 
bitals, values range from 0.4 eV to 0.7 eV at ambient 



9 




K rMKH ALH K TMKH ALH K TMKH ALH 

(a) 96.5GPa - B8 (b) 120GPa - B8 (c) 160GPa - B8 



FIG. 4: (color online) LSDA band-structure of paramagnetic MnO in the B8 high pressure phase. Shown are the Kohn-Sham 
bands (solid) and the eigenvalues of the Wannier Hamiltonian (dashed). The vertical axis is energy (in eV) with the Fermi 
level at the origin. 



B. MnO at high pressure : the paramagnetic B8 
phase 

1. The band- structure 

The B8 phase of MnO is of NiAs (B8) structure. As 
discussed in Ref.^^, the B8 and dBl structure are rather 
similar, and a transition from the high pressure dBl to 
the B8 phase requires a moving of oxygen atoms from 
"distorted octahedral" to "perfect trigonal" positions. 

Yet, in the pressure dependence, there is an impor- 
tant difference : Augmenting pressure in the dBl phase 
leads to an increase of the trigonal angle, and thus to 
anisotropic changes. This results in an almost pressure- 
insensitive bandwidths, as discussed above. In the B8 
phase, the c/a ratio basically remains unchanged, and it 
is thus only the uniform modification of the lattice con- 
stant that accounts for the change in volume. We note 
that the different coordination causes the d-orbitals to 
split into tvifo doubly degenerate e^'^ and a single ag or- 
bitals (see e.g. Ref.^^). Fig. 4 shows the band-structure 
of B8 MnO for various pressures, along with the bands in 
the Wannier gauge for the d-model. The corresponding 
DOS is displayed in Fig. 5. The lowest indicated pres- 
sure is in fact smaller than the largest value we had cho- 
sen for the dBl phase. This is justified by the fact that 
both phases are coexisting over a finite pressure range, 
as evidenced from optics'^'''^^ and x-ray experiments.^^ 

Of course, the LDA yields a metal for all pres- 
sures. Since the Mott transition within the paramag- 
netic B8 phase only takes place at finite temperatures, 
this is congruent with experiments. However, we used 
again experimental structures and volumes at ambient 
temperature,^^ since our aim is the construction of a 



model for the metal-insulator transition within the B8 
phase. Continuing our remarks from the Methods sec- 
tion, this model construction is approximate. While it 
is perfectly admissible that the one-particle dispersion is 
metallic, one might object that the pd screening that we 
include in the determination of the Hubbard U could be 
overestimated since the exact Green's function will have 
excitations with a larger separation because of the insu- 
lating nature of the system. 

In the LDA band-structure. Fig. 4, the d-bandwidth 
becomes slightly larger with increasing pressure. As men- 
tioned in Ref.^^ this restricted increase is a result of the 
fact that the d-bandwidth is mainly constituted by the 
hybridization contribution of the 02p orbitals. Yet, while 
tpd increases with pressure, so does the charge trans- 
fer energy Ac*, and thus the effective hopping, tp^/Act, 
does not increase much (cf. our discussion of the dBl 
phase.). The preponderant effect of pressure on the band- 
structure is a moving away from the Fermi surface of the 
occupied 02p and the unoccupied Mn 4s, 4p and Ad or- 
bitals, owing to greater hybridization splittings. 

2. The interaction 

Given the above changes in bandwidths, the bare 
Coulomb interaction is expected to be roughly pressure 
independent. As can be inferred from Fig. 6 (a), which 
displays this quantity for both the pd and the d-only 
model, this is found indeed to be the case. The absolute 
values are comparable in size to those of the dBl phase. 

The pressure dependence is decisively different for the 
partially screened interaction -the Hubbard U- shown 
in Fig. 6, as above, again for the three different types 
of models. As pressure augments, the changes induced 
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FIG. 5: (Color online) LDA density of states (DOS) of the 
paramagnetic B8 phase of MnO for various pressures. The 
Fermi level corresponds to a; = 0. 



to the bands-structure lead to a decrease in screening 
strength, resulting in an almost isotropic increase of the 
effective interaction of up to about 0.1 eV/lOGPa for the 
pd-model. 

The average nearest neighbor interaction elements Upd 



and (for definitions see above) for the pd-model are 
compiled in Tab. IV. The trend with pressure is the same 
as for the intra-atomic elements, i.e. they grow with pres- 
sure. However, like in the dBl phase, they are smaller 
by a factor of about four with respect to the intra-atomic 
elements. While the pressure dependence of the diago- 
nal terms of the screened interaction could be attributed 
almost exclusively to changes in the screening, it is no- 
table that the off-diagonal elements of already the bare 
Coulomb matrix elements change by almost 10% from 
the lowest to the highest pressure considered. 

What had been said about the Hund's coupling of the 
dBl phase is also true here, as seen from Table VI. Being 
less susceptible to screening processes, it does not show 
a significant pressure dependence. However, it is, again, 
sensitive to the combination of orbitals involved, and val- 
ues can vary by as much as 65%. Since the physics of the 
high spin to low spin transition^^ within the B8 phase is 
sensitive not to the diagonal elements of the interactions 
(which change considerably), but mostly to changes in 
J^^ (that are small), arguments made about the transi- 
tion mechanism^^ may not be affected when using the 
pressure dependent matrix elements instead of constant 
ones. Yet, spectra will change on a quantitative level. 



96 (160) GPa 


dxy 






dyz 








0.49 (0.51) 


0.70 (0.69) 


0.61 (0.61) 


0.58 (0.54) 


dxz 


0.49 (0.51) 




0.48 (0.46) 


0.63 (0.65) 


0.62 (0.63) 


d-,2 


0.70 (0.69) 


0.48 (0.46) 




0.48 (0.49) 


0.69 (0.69) 




0.61 (0.61) 


0.63 (0.65) 


0.48 (0.49) 




0.49 (0.51) 




0.58 (0.54) 


0.62 (0.63) 


0.69 (0.69) 


0.49 (0.51) 





TABLE VI: Hund's coupling J = H^ro o ('^ = 0, a / /J) of paramagnetic B8 MnO for the d-only model and the pressures 96 
and 160GPa. 



IV. SUMMARY, CONCLUSIONS AND 
PERSPECTIVES 

In conclusion, we have studied the impact of exter- 
nal pressure onto the construction of effective low en- 
ergy many-body models for the realistic example of man- 
ganese monoxide. This compound is particularly chal- 
lenging for a first principle description, since energies, 
such as the transfer integrals, the Hubbard U and the 
charge transfer energies are of the same magnitude. On 
the level of band-structures, we concluded that the pres- 
sure induced vanishing of the gap in the antiferromag- 
netic phase of MnO is driven by the crystal field rather 
than by an increase of band-widths, the latter of which 
is found to be limited by the growing rhombohedral dis- 
tortion. The high pressure paramagnetic phase, on the 



contrary, exhibits the usual pressure caused increase of 
hybridizations and bandwidths. Using the constraint 
RPA approach within setups of maximally localized Wan- 
nier functions, we further found that, for a quantita- 
tive description, the pressure dependence of the Hub- 
bard U can not be neglected. While the pressure sen- 
sitivity in the low pressure antiferromagnetic phase of 
MnO is dominated by compression induced changes of 
the bare Coulomb interaction, the effective interaction 
of the high pressure paramagnetic phase is mostly deter- 
mined by modifications in the band-structure that affect 
the screening properties of the system. Our ab initio 
constructed low energy models can provide a quantita- 
tive starting point for the application of sophisticated 
many-body methods. 
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FIG. 6: (Color online) local Coulomb interaction of the paramagnetic B8 phase of MnO for different pressures, (a) diagonal 
elements of the local, bare interaction VJ^q'™™ for the pd and d-only model, and resolved for the different orbitals a. (b) zero 
frequency limit of the RPA partially screened local interaction Uaa = Wro'o ""('^ = 0) for the pd model (with and without p 
screening), and the d-only model. 
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energy system - a reason why the constraint LDA tech- 
nique often overestimates the values of U. 
This discards correlation induced changes in the charge 
density. See e.g.®* for a charge self-consistency scheme 
within DMFT. 

Still, both LDA and RPA work at T = 0. 
In case the bands of the "low" and "high" subspace are 
entangled, the separation P = Piow + Pr is not uniquely 
defined. For a way to overcome this arbitrariness see the 
recent®^. Since in our case the entanglement is not too 
pronounced we do not employ this refinement. 
While the RPA screening and thus also the effective in- 
teractions are energy dependent, so far, most many-body 
techniques only use static interactions, which is why we 
will discuss in this work only the low energy limit, i.e. for 
example the Hubbard Uap = W^r^.o'''^ = 0) (using the 
conventions of Eq. (2)). For a discussion on the frequency 
dependence of the screened interaction see the Ref.^. 
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mentioned combination of GW*®"** with the dynamical 
mean field theory (DMFT)^^ "GW-^DMFT" in which 
the Hubbard U is promoted from a parameter to a quantity 
that obeys a self-consistency relation. 

In this definition, it also includes the respective contribu- 
tions from the exchange splitting. 
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